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Abstract 


THE  BOUND  ARY-INTEGRAL  METHOD 
FOR  PLANAR  MICROSTRIP  CIRCUITS 


X 

The  main  subject  of  this  thesis  is  to  develop  a  general  purpose  computer 
program  for  the  analysis  containing  a  planar  microstrip  circuit  arbitrarily  shaped. 
Based  on  the  boundary-integral  method  (also  called  the  contour-integral  method),  the 
analyses  of  circuits  containing  right  -  angled  microstrip  bend  with  or  without  a  miter 
cut  and  circuits  of  a  T-junctions  with  or  without  a  V-shaped  cut  are  presented  in  this 
thesis.  Scattering  parameters  and  other  related  data  of  the  circuits  are  calculated  and 
presented.  They  agree  well  with  published  data.  It  is  shown  that  the  boundary- 
integral  method  is  an  applicable  and  effective  computer-aided  design  tool  for 
microstrip  circt  rs  of  ..non  complicated  geometries. 
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Introduction 


In  the  development  of  microwave  integrated  circuits  (  MIC )  and  millimeter- 
wave  integrated  circuits  ( MMIC ),  the  planar  microstrip  line  technique  gained  much 
attention  mainly  due  to  the  fact  that  the  mode  of  propagation  on  the  microstrip  line  is 
almost  TEM.  The  so  called  waveguide  model  has  been  applied  successfully  to  a 
number  of  microstrip  discontinuities.  Microstrip  circuits  are  often  accompanied  by 
discontinuities  of  one  type  or  another,  of  which,  the  most  common  forms  are  the 
right-angled  bend  and  T-junction.  The  geometries  of  the  circuits  analyzed  are 
shown  in  Figures  1  and  Figures  2.  In  order  to  simplify  the  analysis  and  to  make  use 
of  well  developed  waveguide  techniques,  the  microstrip  line  circuits  are  usually 
replaced  by  equivalent  parallel-plate  waveguides  with  magnetic  side  walls  (as  shown 
in  Fig.3).  After  each  microstrip  line  network  is  replaced  with  a  closed  structure,  the 
microstrip  discontinuity  problem  becomes  that  of  a  closed  waveguide.  Due  to  the 
discrete  nature  of  modal  spectrum  in  closed  structure,  the  analysis  of  the  scattering 
from  the  discontinuity  is  handled  rather  easily,  typically  by  a  mode  -  matching 
method. 

There  are  many  methods  that  can  be  used  to  analyze  planar  circuits,  such  as 
eigenvalue  method,  or  the  Rayligh-Ritz  method.  For  planar  circuits  with  simple 
geometrical  shapes,  like  circles  or  rectangles  shapes,  the  eigenvalue  method  can  be 
readily  used. 
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The  Rayleigh-Ritz  methcxi  or  finite  element  method  on  the  other  hand  can 
be  used  to  solve  more  complicated  planar  circuits.  However,  the  computation  of 
eigenfunctions  over  the  entire  area  of  the  circuit  pattern  is  time-consuming.  The 
segmentation  and  desegmentation  methods  [  1][  2]  can  be  used  to  analyze  iiregular 
complicated  circuits,  but  there  are  some  limitations  on  the  shapes  of  such  circuits. 
Details  of  these  methods  of  analysis  as  well  as  their  limitations  were  discussed  in 
more  detail  in  a  number  of  articles  [1],[2],[3],[  4]. 

For  an  arbitrarily  shaped  planar  circuit,  using  the  boundary-integral  equation 
method  to  develop  a  general  purpose  computer  program  to  do  the  analysis  is 
effective,  accurate  and  adaptable  since  only  the  voltage  and  the  current  along  the 
periphery  of  the  circuit  --  not  the  entire  area  of  the  circuit  -  need  to  be  considered. 
Circuit  parameters  can  be  derived  from  the  boundary-integral  equation  method 
directly.  Hence,  analysis  of  more  complicated  circuits  can  be  handled  easily  and  the 
computer  rime  is  reduced 

In  the  next  section,  this  boundary-integral  equation  method  is  discribed.  The 
process  of  program  design  is  also  given.  In  the  subsequent  section,  results  based  on 
the  formulation  are  presented  in  terms  of  the  scattering  parameters,  and  a  comparison 
of  results  obtained  by  different  numbers  of  sampling  points  along  the  periphery  is 
also  presented. 
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Formulation 


Based  on  the  two-dimensional  planar  circuit  presented  and  introduced  in  [5],  a 

planar  microstrip  model  with  an  arbitrary  shape  is  considered.  It  is  filled  with  material 
of  effective  permittivity  and  the  effective  width  is  .  There  are  several 


coupling  ports  located  around  the  planar  circuits.  We  divide  the  circuit  periphery  into 

N  incremental  sections  numbered  as  1,2,3 . N,  having  width 

Wi,W2,W3, . ,Wn,  respectively.  Coupling  ports  are  assumed  to  occupy  two  or 

more  sections  (as  illustrated  in  Fig.4).  Sampling  points  are  set  at  the  center  of  each 
section,  and  the  magnetic  and  electric  field  intensities  arc  assumed  to  be  constant  over 
each  section. 

Using  Weber's  solution  for  cylindrical  waves[6],  the  potential  at  each 
sampling  point  along  the  periphery  is  satisfied  by  the  following  matrix  equations  : 


N  N 

luy  Vj  =  XhijIj 

j=1  ^  ■’  j=l 
Mj  =  l  (i=j) 

=  ^  !  cos9yH5^\kr)  ds 


(1) 


(2) 


o)fxd  1 


0)^d 

4 


kW: 

[  I-2i-(ln-^-l+y)] 

71  ^ 


(i=j) 


(3) 
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Y  =  0.5772 . (  Euler's  constant )  k=0) 

where  I  =  -2^  W. ;  the  total  current  flowing  into  the  jth  port  on  both  the 

upper  and  lower  surfaces  of  the  circuit  plates. 

/O)  (0) 

.•  the  zeroth  order  and  the  first  order  Hankel  function  of 

the  second  kind,  respectinely. 

d :  the  thickness  of  dielectric  layer . 
r  and  :  the  distance  between  the  ith  and  jth  coupling  ports  and 

the  angle  made  by  the  line  connecting  these  points  and 
the  normal  at  the  jth  sampling  point . 


The  rf  voltage  of  each  sampling  point  can  be  obtained  by  solving  (1): 


[V|  =  U'‘h(I] 

(4) 

[V]  =  [ZU  [1 1 

(5) 

5 


^2 


Vn 

\  / 


/ 


,11 

21 


nl 


[Zl.xn  =  U'*H 


^In 
Z 

2n 


\  /  \ 

I 


nn 


I 


1 


\  V 


(6) 

(7) 


[Z]  is  an  nxn  matrix.  It  is  considered  that  all  n  sections  are  coupling 

ports.  In  fact,  most  of  the  n  ports  described  above  are  open  circuits;  in  other  words, 
most  of  the  nxn  elements  of  are  zero.  Therefore[Z^^j^,  can  be  reduced  by 

choosing  only  the  necessary  elements  given  by  (7).  The  reduced  matrix  is 

described  in  (8).  N  is  the  sum  of  all  sampling  points  along  all  A  ports.  (As  shown  in 
Fig.5) 


A 

N=X 

i=l 

m  :  the  number  of  sampling  points  of  ith  port 


[V]=[ZkN[I] 


(8) 
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(9) 


The  next  step  is  to  consider  the  higher  order  modes  of  each  port  by 
expanding  the  RF  voltage  and  current  density  of  each  port  in  terms  of  m  stripline 
modes. 


=  I  ^  L  COs[(L-l)^] 

L=1  Wj 

I(xjib=  Z  - '0{^^COs[(L-l)-^] 

L=1  m;  W. 


(10) 


(11) 


:  kth  sampling  point  at  the  ith  port 
k  =  1,2, . ,n^  i  =  1,2, . ,A 

%  magnitude  of  modal  voltageat  each  sampling  point 


rr^ :  the  total  number  of  sampling  points  at  the  ith  port 

(i  =  1,2 . ,A) 

Wj  :  width  of  the  ith  port 


From  Fig.5,  we  obtain  the  point 

Wj  Wj(k-l) 


x(i)_ 


2nij 


m: 


w- (2k-l) 

2mj 
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So,  equations  (10)  and  (1 1)  can  be  witten  as  follows: 


V(x®)=  ] 

A  ^  201; 


I  (X®)  =  I  — a>{;>cos[  ] 


L=1  m 


and  then 


V(X^/>)  =  — )+. 

2m^ 


(1)  (mj-l)7l 
.+«  cos(— - ) 


v(xs,")=u<,‘>*Dy'cos(2i:ii^)*. 

2m<i 


(1)  (2mj-l)7l 

V(X^  )='^i  ^2  ^cos( - )+. 

*1  7m. 


(j)  (2k-l)(mj-l)7l 

.+x>  jj,  cos( - ) 

1  2mi 


(2nij-l)(mj-l)7C 
.+\)  „  cos( - ) 


/AY  /AY  /AY  K  (A)  (ni.-l)7C 

V(X<'‘>)  =  D<''>*uWcos( - h . +1)  cos(— ^ - ) 


v<W»s(^^). . 


(A)  (2mA-l)7C 

V(Xm  )=d;aIu^a)cos( - )+. 


(A)  (2mA-l)(niA-l)7C 
.+1)  ^  cos( - 


KX^)  =  I  i^2  ^C0S(  — )+. 

m,  2ni, 


i(x<;))=J_(ia),ia)cos(SiH)^, 


<’>  (1)  ,n  (2m|.i)7t 

>='i  ^‘'cos( - )*. 


(1)  (mj-l)7C 


(1)  (2k-l)(lTl|-l)n 

.+i  COS( - ) 

1  2m, 


(2mj-l)(mj-l)7l 
.+i  m  cos( - ; 

I 


i(A).2)cos(  -^)+. 
niA  2mA 

I(X<f  >)=-L  (i(A),i(A),^3(£l^)^ 


(A)  , 


(2mA-l)7C 


(A)  (m^-l)7C 

.+i  m  cos( - )\ 


(A)  (2k-l)(mA-l)7C 

.+ii^  cos( - )) 


(2mA-l)(mA-l)7l 
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(Nxl) 


(Nxl) 


12 


’  . l5r> 


*  . cos(3(ni-I)  S— ) 


1  cos((2m-l)  ^ — ) 
1 


1  '  2in, 


cos((2m-lXm-l) 

1  1  -inii 


^  . 

‘  “"5  is-'  . cosp(m,-l)  Ij.) 


^  cos((2m-l)  -S — ) 

i  2mi 


cos((2m  -lXm  -1) 
1  1 


*  . 

'  . cosp(m-l)  1^) 


^  cos((2ni  -1)  5 — ) 
A  2iny^ 


cos((2m  -IXm  -1)  -5—) 

A  A 


(1) 

M  = 


m  X  m 
1  1 


We  define 


m  =tzj^[i] 

[A)[-U)  =  [Z]  [A][M)[i] 


[D  ]  and  [  i  ]  are  Nxl  matrices 
[  Z  ],  [  A  ]  and  [  M  ]  are  NxN  matrices 


[B]  =  [Z][A][M] 


[A][\)]=  [B]  [i] 
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For  the  purpose  of  separating  the  first  order  mode  from  higher  order  modes, 
the  matrices  [  A  ],  [  B  ]  are  rearranged  as  below: 


5 


B  B 

IJ'f-m  +1  12 

!  A 


B  i 

N-n^+i‘i‘ 


A 


N-n^+lJ^-m^+1  N-nj^+1,2 


B 

N-nj^+lJ4 


B21  B2,m+l**“®2J^-m+l  ®22  ••••®2m‘ .  ®2 +2  •  ®2N 

:  1  .  A  .  .1  A 


m  1  ®m  J>I-in  +1  ®in  2""®in  m 

:  1  ;  1  A  1:11 


®N-m+2,l  ®N-in +2.N-m +1  ®N-m +23? 
A  -A  A  .A 


®N1  . ?NJNr-m  +1  ®N2 

A 


®N-m +2j^-ni*+2  ®N-m  +24v; 
.  A  A  .  A 

Bvt  kt  _  .  ^  BvrikT 


• 

(1) 

1 

2 

■(1) 

i 

mi 

(A) 

i 

7 

*Na^-m  +2 
A 


0 


^12. 


0 


(2) 

11 


0 


^Im. 


,{2) 

^12.. 


A  (2) 

•  ^Im. 


-(A) 

^12.. 


A(A) 

^Im. 


''aS*  0..: 

\ 

^22  .. 

^2m 

0 . 

\ 

..0 

aSo... 

0  ’a® 

0 . 

^nj2  . 

••  TlAg 

aI^) 

••  ^.^2 

0 

C21  = 

o'  a(2^ 

0 . 

^-22 

0 . 

aO.) 

\2 

X(2) 

-  Tz 

0.' . 

oaJJ) 

0 . 

A  (A) 
^22  .. 

^2m, 

A 

0 . 

o 

> 

3 

>Z^ 

0 . 

< 

a(A) 

A^2. 

A'CA) 

•  •n/ 

/b  b 

I  ,11  l^^+l 


B  ^ 

i;^-m  +1 
*!  A 


B  - 

\^N-nj^+l,l 


B 

N-m  +1  Jl-m  +1 J 
A  A  X' 
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An  important  property  of  waveguide  modes  is  the  existence  of  a  characteristic 
wave  impedance.  From  the  field  theory  point  of  view,  each  possible  solution  to  the 
wave  equation  is  called  a  mode.  A  TEM  wave  is  called  a  transmission-line  mode, 
and  all  other  waves  are  called  higher-order  modes. 


19 


When  assuming  that  the  stripline  is  straight  with  sufficient  length,  only  the 
TEM  mode  can  propagate  in  the  uniform  stripline,  and  higher-order  modes  which 
exist  at  the  port  i  are  considered  to  be  loaded  by  the  reactive  characteristic  impedance 
matrix  [ 

L  =  2,3 . m^ 

i  =  1,2, . ,A 

-jZj((L-l)X/2W^-y^)^-l)-l/2  (p  =  q)  (26) 

=  0  (27) 

p,  q  =  1,2,3, . ,m  -1  i  =  1,2, . ,A 


1207td 

'\j^c 


d  ;  the  thickness  of  the  substructure 
y/'  the  effective  width  of  the  stripline 

:  the  effective  relative  dielectric  constant 
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(N-A)x(N-A) 


N  =  m|+n^+ . 


(29) 


(30) 


Then,  we  can  find  the  relative  equation  for  the  higher  order  modes  at  each  port. 


V2=2lI2 


(31) 


Based  on  analysis  using  contour-integral  method,  (shown  as  Eq.44) 

the  fundmental  modes  of  the  network  at  pon  1,2,....,A  can  be  extracted  from  TEM 
modes  and  loaded  by  the  characteristic  impedances  of  the  higher  order  modes.  [Z]^^ 

can  be  transformed  into  the  scattering  matrix  ( shown  as  Eq.46)  directly  as  follows  . 


CnVi+q2V2= 

C21V,  +^2^2=  P21I1  +0.2^2 
^2=  Vl 


(32) 

(33) 

(34) 


^11^1  q2^^2  ~  ^2^ 


(35) 
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QlVj  =  E^llj  +  (  1^2’ ^2^)  4 

c„v, +(qj-D,2ZL')  v^  =  1:1,1,  (39) 

C21V,  +  ( q^ .  q^zL')  Vj  =  (40) 

from  (39), (40) 

(  q^.  q^zL')  Vj=  - QiV,+  q.i,  (4i) 

Vj=.(qj.qjZL')''c2iV,»(  q^. q^zLY’D^,!,  (42) 

substitute  (42)  into  (39) 

[ -  ( c;2  -  DijZi.')  (  q2  -  q2ZL)'*C2,]  v,  = 

[qi-(  qj-qjZL'x  q^-q^zLY’Di,]],  (43) 
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Numerical  results 

1.  Right-angled  bend  with  and  without  a  miter  cut :  ( Fig.l  a.b ) 

Shown  in  Fig.6  is  the  variation  in  the  reflection  coefficient  I  for 

chamfered  and  unchamfered  right-angle  bends.  The  results  are  compared  with 
Gupta's  results  [9].  In  order  to  reduce  the  effect  of  the  discontinuity  reactances,  an 
isosceles  triangle  part  must  be  removed.  The  size  of  the  triangle  removed  has  been 
optimized,  to  minimize  the  reflection  coefficient  .  Figure  7  shows  that  when  a 
measured  cut  ratio  c/w  is  equal  to  about  0.8,  the  V.S.W.R.  is  minimized.  Kaneko 
has  the  same  conclusions[7]. 

2.  T-junction  without  a  miter  cut :  ( Fig.  2a ) 

Results  for  frequency  dependent  scattering  matrix  coefficients  for  a  T- 
junction  (^=  9,7,  d  =  0,0635  cm,  v^=0.056  cm )  are  shown  in  Fig.8.  It  notes  that 

I  S22I.  as  shown  in  Fig  8a,  the  power  reflected  by  the  T-junction  discontinuities 
increases  with  frequency  in  the  range  0<  f  <  f  j,  where  is  the  cut-off  frequency 

of  the  first  higher  order  mode.  Correspondingly,  the  power  transmitted  from  port  1 
to  port  2  decreases  with  frequency  greater  than  f  j.  The  power  of  a  port  is  transmitted 

by  the  first  higher  order  mode,  so  that  the  transmission  coefficient  I  S|2l  is  alw?ys 
smaller  than  that  for  f  =  0  ( shown  in  Fig.8b ). 
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Fig.Sc  indicates  that,  if  the  TEM  mode  is  incident  at  port  1,  the  reflected 

power  decreases  with  increasing  frequency,  and  another  power  transmitted  to  port  3 
behaves  in  an  opposite  sense  to  I  Si2l'  Also  plotted  in  Fig.8  are  the  numerical 

results  reported  in  [8]. 


3.  T  -  junction  with  a  V  -  shaped  miter  cut :  ( Fig.  2b ) 


This  section  considers  a  T-junction  with  an  isosceles  triangle  with  angle  6  = 
45°removed;  the  side  of  the  triangle  removed  is  equal  to  c  (cut  depth).  For  a  cut  ratio 

c/w2  =  0.4035...,  the  scattering  coefficients  are  plotted  in  Fig.9.  The  reported  data  [9] 

is  also  plotted  in  fig.9.  When  the  frequency  increases,  the  variation  of  the  magnitude 
of  reflected  coefficient  l^jl  keeps  almost  constant  and  is  between  a  range  of  0.41  to 

0.45.  The  predicted  value  of  1  %2l  is  greater  than  the  published  value  by  2.5 
percent. 

For  different  cut  ratios.  Fig.  10  shows  clearly  that  by  removing  an 

isosceles  triangle  from  a  T  -  junction  with  0  =  45Ss(o, )  and  c/w2  =  0.8,  the  effects  of 
I  Sjjl  are  minimized  at  the  higher  frequency  band .  This  phenomenon  is  similar  to  the 

case  of  right-angled  microstrip  bend  with  a  miter  cut  ratio  0.8.  By  considering  a  plane 

of  symmetry  located  at  half  the  width  of  port  2,  the  T-junction  can  be  compared  with 

a  right-angled  bend  on  just  one  side  of  the  symmetrical  plane.  The  reflection 
coefficient  I  is  equal  to  I  ^3!  of  the  corresponding  side  of  the  bend.  Thus,  if  only 

I  S|  j  I  or  I  ^3!  is  to  bt  optimized,  then  only  a  right-angled  bend  with  one  half  of  the 


width  of  port  2  need  to  be  considered. 
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In  the  boundary  -  integral  method,  the  circuit  periphery  is  divided  into  N 
sections  ,  the  value  of  N  determines  the  size  of  NxN  impedance  matrix  [Z]  .  The 
greater  N  is  ,  the  more  precise  the  result  will  be.  The  error  decreases  from  8  to  2.5 
percent  with  N  increased  from  80  to  300  ( Fig.l  1).  Also,  for  the  purpose  of  saving 
computer  time,  the  problem  of  properly  selecting  the  location  of  each  scattering  point 
become  more  important  when  the  effective  width  of  each  port  is  not  the  same. 
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Conclusion 

Computer-aided-design  is  a  very  powerful  tool  in  the  circuit  design  of  MIC 
and  MMIC.  A  great  savings  of  time  and  money  can  be  resulted  if  a  design  is  verified 
before  it  is  fabricated.  The  work  described  in  this  paper  showed  how  the  boundary- 
integral  method  can  be  used  to  solve  planar  microstrip  circuits.  The  results  discussed 
here  are  in  excellent  agreement  with  other  methods  for  the  right-angled  bend  and  T- 
junction.  Future  studies  will  concentrate  on  developing  this  method  for  more 
complicated  microstrip  discontinuities. 
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